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SUMMARY 


An  experiment  has  been  performed  which  demonstrated  that  electron 
temperatures  and  plasma  effects  similar  to  those  predicted  for  electro- 
magnetic pulses  (EMP)  propagating  in  the  high  altitude  source  region  and 
D-region  of  the  ionosphere  can  be  maintained  and  observed  in  a laboratory 
experiment.  Scaling  laws  have  been  developed  which  allow  a long  propaga- 
tion path  to  be  modeled  by  a one  meter  TEM  waveguide  emersed  in  a 
laboratory  plasma.  It  has  been  demonstrated  that  nonlinear  collision 
absorption  processes  and  electron  avalanche  do  occur  as  qualitatively 
predicted.  The  data  collected  in  this  experiment  shows  that  absorption 
of  energy  from  propagating  signals  can  occur  and  frequencies  appreciably 
above  the  plasma  frequency. 
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PREFACE 


The  theoretical  work  underlying  the  experimental  effort 


described  in  this  report  has  been  developed  by  Capt.  William  A. 


Seider  at  the  Air  Force  Weapons  Laboratory,  Kirtland  AFB,  New  Mexico. 


The  plasma  machine  and  the  experimental  program  herein  described  has 


proceeded  under  the  direction  of  Dr.  Robert  N.  Carlile,  Professor  of 


Electrical  Engineering,  Department  of  Electrical  Engineering,  Uni- 


versity of  Arizona,  Tucson,  Arizona.  Three  graduate  students  have 


been  involved  with  the  development  of  this  equipment:  Robert  A. 


Piejak,  PhD  candidate;  Gary  L.  Jackson,  PhD  candidate,  and  Gary  L. 


Hagedon,  MS  candidate. 
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I.  Introduction 


High  and  low  amplitude  electromagnetic  pulse  (EMP)  signals  can 
propagate  into  space  and  pose  a potential  threat  to  satellites.  At  the  Air 
Force  Weapons  Laboratory, (AFWL) , a theory  has  been  developed  which  predicts 
that  a large  amplitude  EMP  signal  propagating  through  the  D-region  of  the 
ionosphere  will  cause  rapid  heating  of  the  electron  in  this  region,  result- 
ing in  conditions  being  such  that  considerable  energy  is  absorbed  from  the 
pulse.  This  absorption  would  thus  reduce  the  potential  threat  of  the  EMP 
to  satellites. 

In  June,  1974,  the  Air  Force  contracted  with  The  University  of 
Arizona  to  build  an  experimental  machine  that  would  allow  the  AFWL  theory 
to  be  tested  experimentally.  During  the  6-month  technical  phase  of  the 
contract,  the  emphasis  has  been,  of  necessity,  on  the  construction  of  the 

machine.  On  January  15,  1975,  when  this  phase  ended,  the  machine  had  been 

\ 

constructed,  and  some  preliminary  data  had  been  obtained  from  it.  This 
report  describes  the  machine  in  detail  and  presents  these  preliminary  data. 

In  Section  II  of  this  report , a summary  of  the  theory  developed 
at  AFWL  is  presented.  In  Section  III,  a general  description  oi  the  experi- 
mental procedure  that  will  lead  to  a testing  of  the  theory  is  given, 
including  a discussion  of  scaling  between  the  ionospheric  plasma  and  the 
laboratory  plasma.  In  Section  IV  the  features  of  the  machine  are  discussed 
in  somewhat  general  terms.  In  Section  V,  the  details  of  the  construction 
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A theory  developed  at  the  Air  Force  Weapons  Laboratory  predicts 
electron  heating  of  the  D-region  by  an  EMP.  A summary  of  that  theory  is 
presented  here^.  Note  that  all  figures  in  this  Section  are  from  Reference 


To  summarize  the  theory  in  words,  consider  a large  amplitude  fast 

risetime  EMP,  incident  upon  the  lower  edge  of  the  D-region  (low  altitude 

edge) . As  the  pulse  starts  to  propagate  through  the  D-region,  it  sets  up  a 

conduction  current  density  J.  The  current  density  J results  in  joule  heat- 

— — 3 

ing  of  the  electrons  at  a rate  of  J • E joules  per  second  per  m . This 
energy,  absorbed  by  the  electrons  from  the  EMP,  is  ultimately  transferred 
to  the  neutrals. 

Joule  heating  of  the  electrons  roughly  increases  with  altitude  in 
the  D-region,  so  that  the  energy  in  the  pulse  declines  at  first  slowly,  and 
then  rapidly  as  altitude  increases.  The  rate  at  which  energy  is  transferred 
to  neutrals  versus  the  rate  at  which  it  is  absorbed  from  the  pulse  is  such 
that  the  electron  temperature  reaches  a peak  around  85  km  which  may  be  of 
the  order  of  several  hundred  electron  volts  if  the  pulse  amplitude  is  high 
enough.  These  high  electron  energies  cause  electron  avalanching  to  occur 
so  that  the  electron  number  density  will  increase.  The  avalanche-produced 
electrons  will  contribute  to  J.  Nearly  total  absorption  of  the  pulse  can 
occur  by  the  time  the  pulse  emerges  from  the  high  altitude  side  of  the  D- 
region. 


The  starting  point  for  the  theory  is  the  first  three  moments  of 
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the  Boltzmann  equation  and  Maxwell’s  equations: 
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Maxwell’s  equations: 


(Third  Moment) 
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The  system  is  shown  in  Fig.  1.  In  Equations  (1)  through  (8),  ng  is  the 

number  density,  v is  the  average  velocity,  E is  the  amplitude  of  the  EMP. 

J = - n e v is  the  current  density,  T is  the  electron  energy  expressed  in 
e e 
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electron  volts,  and  B is  the  magnetic  field.  These  quantities  are  all 
functions  of  the  retarded  time  t only.  B can  contain  a DC  part  due  to  the 
earth's  magnetic  field.  All  quantities  except  are  expressed  in  terms  of 
their  usual  MKS  units.  When  v^,  v^,  and  are  multiplied  by  o,  the  number 
density  of  neutrals,  the  results  are  collision  frequencies.  These  are 
associated  with  collision  cross  sections  for  electron  generation  through 
avalanching,  momentum  transfer,  and  energy  transfer  to  neutrals,  respectively. 
They  are  functions  of  the  electron  temperature.  Tgo  is  the  ambient  electron 
temperature.  Other  constants  should  be  self-evident. 

The  transformation  to  retarded  time  is  made, 

tg  = t - r/c  (9) 


r 21 

where  t is  clock  time  and  c is  the  velocity  of  light  . 
transformation,  3/3t  is  replaced  by 


According  to  this 


and  V is  replaced  by 


V - r (1/c)  3/3t 


where  r points  from  the  source  point  to  (r,9,$).  [See  Fig.  1.] 

After  this  transformation  has  been  made,  it  is  possible  to  show 
that  variation  of  the  dynamical  variables  with  r is  much  slower  than 
variation  with  cts-  We  can  then  neglect  spatial  derivatives  compared  to 
derivatives  with  respect  to  ts-  Equations  (1)  through  (8)  have  been  trans- 
formed according  to  the  procedure  described  above. 

The  D-region  is  subdivided  into  incremental  lengths  along  the 


r axis  in  Fig.  1 of  10  km  each,  corresponding  to  an  incremental  altitude 
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change  of  about  1 km.  An  ambient  number  density  and  temperature  is 
assigned  to  each  increment.  An  electric  field  (Fig.  2)  is  applied  to  the 
input  to  the  first  increment,  which  occurs  at  52  km.  Four  different  cases 
have  been  considered  at  AFWL,  indicated  by  the  four  pulses  labeled  A,  B, 

C,  and  D in  Fig.  2. 

A computer  iterative  technique  was  used  to  determine  the  response 

of  the  plasma  to  the  pulse,  and  effect  of  the  plasma  to  the  pulse.  Consider 

an  increment  whose  input  is  at  r^  and  output  is  at  r^^.  Call  this  the  nC^ 

increment.  By  keeping  the  electric  field  at  the  input  to  the  previous  six 

increments,  a Taylor  series  could  be  constructed  to  determine  the  electric 

field  at  r . This  electric  field  is  then  used  in  the  plasma  Equations  (1) 
irri 

through  (3)  as  a source  to  find  v,  T and  a as  a function  of  t . The 

e e s 

boundary  conditions  on  these  equations  are: 


At 


t.-°» 


v = 0, 


T = T , 
e eo 


n = n 
e eo 


From  this,  J is  computed,  and  using  (4)  through  (6),  a correction  to  E at 

r^^  is  found.  This  process  is  repeated  using  the  corrected  E such  that  a 

further  correction  in  E is  calculated.  This  procedure  repeats  until  E at 

r , , stabilizes. 
n+1 

The  results  of  this  procedure  applied  to  nitrogen  ionosphere  can 
be  summarized  as  follows.  The  percent  of  total  energy  stored  in  the  pulse 
as  a function  of  altitude  is  shown  in  Fig.  3.  Most  of  it  remains  in  the 
pulse  until  about  80  km,  and  thereafter,  it  is  rapidly  absorbed  by  the 
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electrons.  The  high  amplitude  pulses  (A,  B,  and  C)  retain  their  energy  j 

£ %■  , 
7, 

' 

longer  on  a percentage  basis  than  the  low  amplitude  pulse  D. 

In  the  wake  of  the  pulse,  the  energy  density  in  the  D-region  (located 

F 

B 

both  in  neutrals  and  electrons)  calculated  at  AFWL  is  shown  in  Fig.  4.  It  | 

i 

generally  increases  with  altitude.  This  energy  is  divided  between  the  i 

electrons  and  neutrals.  The  rate  at  which  energy  is  flowing  into  the  ( 

Ijr 

* i 

neutrals  is  shown  in  Fig.  5.  1 

Energy  flows  to  the  electrons  from  the  pulse  at  a greater  rate  than 

Bi!  ’ 

it  is  absorbed  from  the  electrons  by  the  neutrals.  This  causes  the  electrons  j 

S 

to  increase  in  temperature.  The  peak  temperature  achieved  by  the  electrons  ; 

is  shown  in  Fig.  6.  It  is  seen  that  for  the  high  amplitude  pulse,  an  elec- 

1 

sj 

tron  temperature  of  nearly  100  ev  is  achieved. 

i 

III.  General  Description  of  Experiment  j 

The  D-region  of  the  ionosphere  is  characterized  by  low  electron  | 

1 1 

E-xi  i 

3 ! 

densities  (100  electrons  per  cm  for  a sunspot  minimum,  night  time  ionosphere) . J 

m ■ 

4 i 

For  a collisionless  plasma,  frequencies  below  9 x 10  MHz  would  be  reflected.  } 

s 

As  we  have  seen  in  the  preceding  section,  the  theory  developed  at  AFWL  has 

1 

sy 

indicated  that  in  the  D-region,  a nonlinear  absorption  can  occur  which 

i 

affects  frequencies  considerably  higher. 

ij 

iff 

We  have  conceived  of  an  experiment  which  allows  pulse  propagation  * 

s&. 

in  the  ionosphere  to  be  simulated,  and  which  will  permit  theoretical  pre-  f 

life 
1 ??•"•' 

§|  3fc 

dictions  of  Section  II  to  be  tested.  In  brief,  a TEM  waveguide  is  immersed  | 

P WJEft  < 

1 

in  a nitrogen  plasma.  We  apply  a pulse  simulating  an  EMP  at  the  input  and  | 

i Ifllli 

1 Iplll 

observe  the  received  signal  at  the  output.  We  measure  the  real  time  varia-  | 

■ 

tion  of  the  input  and  output  pulses  with  a sampling  scope. 

6 

6 


BSS5E 


:••  5 


r-  ■ I 


« *-} 


The  TEM  waveguide,  which  is  about  1 meter  long,  simulates  an 
incremental  distance  in  the  D-region  of  the  order  of  1 kilometer.  This 
would  correspond  to  one  iteration  interval  of  the  AFWL  theory. 

We  can  measure,  by  means  of  an  electrostatic  analyzer  probe,  the 
electron  distribution  function.  This  allows  us  to  measure  the  maximum 
energy  of  electrons,  which  the  theory  predicts  can  exceed  100  ev.  In 
addition,  in  principle,  we  should  be  able  to  obtain  from  the  distribution 
function,  the  final  values  of  the  electron  density  and  drift  velocity  after 
electron  heating  by  the  EMP. 

Ambient  values  (pre-pulse)  of  electron  density  and  electron  tempera- 
ture are  measured  by  standard  electromagnetic  cavity  perturbation  and  Lang- 
muir probe  techniques. 

We  need  to  discuss  the  correspondence  between  the  behavior  of  the 
pulse  propagating  in  our  laboratory  plasma  and  that  of  a pulse  propagating 
in  the  ionosphere.  We  make  a correspondence  between  our  1 meter  length 
plasma  filled  waveguide  and  an  incremental  distance  in  the  ionosphere  dis- 
cussed in  Section  II.  The  laboratory  plasma  is  uniform  and  we  assume  that 
the  ionospheric  plasma  over  the  distance  increment  is  also  uniform. 

We  now  discuss  the  determining  Equations  (1)  through  (8)  without 
reference  to  a particular  system.  It  is  assumed  that  these  equations  are 
valid  for  either  system. 

The  uniform  plasma  is  characterized  by  an  ambient  number  density 
ng^  <md  electron  temperature  Tgo  electron  volts.  From  these,  we  may  con- 
struct the  following  constants: 

(1)  The  ambient  thermal  velocity: 


V.  3 £e  T /■] 

t eo 


1/2 


(12a) 


m 


(2)  The  ambient  plasma  frequency: 


and 


f 2,  ,1/2 
a)  = In  e /me  I 
p eo  o 


(3)  The  electron  cyclotron  frequence 


to  = e B /m 
c o 


derived  from  the  DC  magnetic  field: 


— A 
B = u B 


(4)  The  debye  length: 


XD  " V“p 


(5)  The  ambient  electron  momentum  collisions  frequency 


m eo 


We  may  now  define  two  normalized  independent  variables 


t = to  t 
p s 


R = r/X 


D 


and  five  normalized  dependent  variables: 


N = n /n 

e eo 


T = T /T 
e eo 


V = v/V 
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which  are  functions  of  t.  We. also  define  the  auxiliary  function 

F,(T  .T)  = v,(T  )/v  (T  ) 
i eo  i e m eo 


and 


F (T  ,T)  = v vT  )/ v (T  ) 
m eo  m e m eo 


F (T  ,T)  = V (T  )/v  (T  ) 
w eo  w e m eo 


The  set  of  Equations  (1)  through  (8)  may  now  be  put  into  normalized 
form  by  multiplying  each  equation  by  the  following  factor: 
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by 
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(7) 

and 

(8)  by 

(e/m)XD/(op 

The  results  are: 
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(Plasma  Equations) 
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The  electric  field  at  the  input  to  the  ionospheric  increment  or  to  the 
waveguide  has  an  amplitude  Eq  and  a rise  time  t and  fall  time  t^.  Suppose 
that  we  require  that  for  both  systems,  the  following  system  parameters  all 


be  the  same : 
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Then  the  normalized  equations  will  look  the  same  for  either  system,  and 
their  solution  will  be  identical  of  either  system.  To  restate:  if  the 

system  parameters  listed  in  (27)  can  be  made  the  same  for  both  systems,  then 
the  two  systems  will  behave  in  exactly  the  same  way.  The  normalized  inde- 
pendent variables  versus  the  independent  variables  t and  R will  be  the  same, 
and  if  this  set  of  normalized  independent  variables  can  be  found,  say,  in 
the  laboratory  system,  the  unnormalized  variables  in  '.he  ionosphere  can  be 
found  from  (14)  through  (18) . 

This  result  is  much  more  general  than  indicated  here.  The  previous 
discussion  has  made  use  exclusively  of  the  retarded  time  and  not  the  usual 
clock  time.  However,  we  have  shown  elsewhere  that  if  one  starts  with  the 
Boltzmann  equation  and  Maxwell ’ s equation  not  transformed  by  retarded  time, 
but  using  usual  clock  time,  the  conclusions  are  the  same;  i.e.,  if  the 
quantities  in  (27)  are  the  same  for  both  systems,  normalized  dependent  vari- 
ables will  behave  in  the  same  way  in  both  systems  versus  the  independent 
variables  t and  R = r/A^.  All  nonlinear  effects  are  encompassed,  as  well 
as  such  phenomena  as  Landau  damping. 

We  consider  next  the  practicality  of  making  the  nine  quantities 

and  functions  in  (27)  the  same  in  the  laboratory  and  ionospheric  systems. 

Consider  the  first  five,  all  containing  u>  , which  tends  to  be  low  in  the 

P 

ionosphere  (90  KHz),  but  which  can  conveniently  be  at  least  9 MHz  in  the 
laboratory.  If  S is  a scale  factor,  i.e.,  the  ratio  of  u>^  in  the  laboratory 
to  that  in  the  ionosphere,  then  in  the  laboratory,  time  tg  is  scaled  down  by 
S including  t^  and  t^  while  Eq  is  scaled  up. 


Specifically,  if  Eq,  t^,  t^,  v,  and  are  the  ionospheric  values, 

the  correct  values  for  the  laboratory  are  SE  , t /S,  t,/S,  Sv,  and  Su>  . In 

o r f c 

the  laboratory,  EQ,tr,  and  t^  can  be  varied  by  design  of  the  pulse  generator 

as  described  in  Section  V,  while  v and  can  be  varied  by  means  of  the 

pressure  and  confining  magnetic  field.  Since  these  quantities  can  be  changed 

over  wide  limits,  the  appropriate  scaled  values  can  be  usually  achieved. 

It  is  the  remaining  four  quantities  which  are  difficult  to  scale. 

V , F , F , and  F depend  on  T . Thus , T must  be  the  same  in  both  systems 
tim  w r eo  eo  

In  our  steady  state  plasma,  T is  about  1 ev  and  cannot  be  easily  controlled 

while  in  the  D-region,  it  is  about  .01  ev.  Thus,  V^/c  will  be  about  a factor 

of  30  larger  in  the  laboratory,  and  the  functions  F , F , and  F.  will  be 

m w x 

somewhat  different  functions  of  T in  the  ionosphere  than  in  -the  laboratory, 
particularly  F..  Thus,  exact  correspondence  of  the  quantities  in  (27)  is  not 
possible  in  the  ionosphere  and  laboratory. 

However,  the  purpose  of  this  work  is  to  test  the  validity  of  the 
theory  in  Section  II.  To  this  end,  the  information  about  scaling  above  can 
be  put  to  the  following  use: 

A computer  run  should  he  made  with  T set  equal  to  1 ev  for  one 

distance  increment  and  a>  ,E  , t , t,,  v,  and  w chosen  such  that  a)  the 

port  c 

scaled  values  are  attainable  in  our  laboratory  machine , b)  the  theoretically 

predicted  results  show  electron  temperature  approaching  100  ev.  The  results 

from  this  run  can  now  easily  be  scaled  to  predict  our  experimental  results 

for  a range  of  u>  in  our  laboratory  plasma . 

P 

In  the  machine  which  we  have'  built,  the  plasma  is  confined  by  a 
magnetic  field  as  described  in  Section  IV.  However,  according  to  the 
Statement  of  Work  for  this  contract,  the  theory  to  be  tested,  initially 


assumes,  zero  magnetic  field.  The  reader  is  referred  to  Item  6.1  and 
Equations  (1)  through  (4)  i :em  3.2  cf  the  Statement  of  Vcrk.  This  pre- 
sents no  problem  at  medium  l high  altitudes  in  the  D-region  since  we  can 

make  u>  large  enough  and  to  small  enough  that  the  ratio  to  /to  in  less  than 
p c c p 

0.1.  Over  the  length  of  the  TEM  waveguide,  such  a low  magnetic  field  will 
have  negligible  effect. 

At  lower  altitudes,  we  could  use  a discharge  plasma  system  which 
has  not  been  built  as  of  January  15,  1975  (but  which  could  be  added  at  small 
expense,  about  $500).  In  this  case,  the  plasma  is  confined  by  a glass  pipe, 
and  no  magnetic  field  is  required.  This  will  be  a truly  zero  magnetic  field 
plasma. 

However,  it  is  quite  easy  to  take  into  account  the  earth's 

magnetic  field.  AFWL  has  completed  a theory  which  includes  this  magnetic 

field ^ . We  must  require  that  the  direction  of  the  magnetic  field  should 

be  parallel  to  the  direction  of  propagation.  We  can  do  this  by  using  the 

*■ 

electron  beam-generated  plasma  system  and  letting  the  magnetic  field  which 

confines  the  plasma  have  an  appropriate  value.  To  use  the  existing  flat 

plate  TEM  waveguide,  the  magnetic  field  must  be  such  that  the  electric  field 

does  not  demonstrate  Faraday  rotation  significantly  over  the  length  of  the 

waveguide.  The  theoretical  data  presently  available  from  AFWL  show  that 

for  typical  magnetic  fields  found  in  the  D-region  this  will  not  happen. 

For  a strictly  zero  DC  magnetic  field,  and  perfectly  conducting 

planes,  the  fields  of  the  pulse  in  the  waveguide  are  exactly  TEM  fields 

with  identical  electric  and  magnetic  field  configurations  and  dispersion 

curve  to  that  which  would  be  found  in  an  infinite  cross-section  of  iono- 
[33 

spheric  plasma  . We  do  not  use  perfectly  conducting  planes,  but  rather 
planes  consisting  of  stainless  steel  wire  mesh.  We  assume  that  the  fineness 


at 

1 - V " - . 
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of  the  mesh  (compared  to  the  shortest  wavelength  component  of  the  pulse) 

P % 

and  its  slightly  lossy  character  do  not  significantly  alter  the  fields  or 

if-  --. 

dispersion  curve. 

v - 

The  presence  of  a DC  magnetic  field  will  cause  the  fields  to 

I ■JS 

rotate  (Faraday  rotation) . For  the  choice  of  parameters  to  be  used  in 

i _ . 

il 

these  experiments  which  are  to  simulate  the  ionosphere,  Faraday  rotation 

' 

is  evidently  very  slight,  as  has  been  demonstrated  by  the  AFWL  theory. 

! 

e.g.,  if  at  the  input  to  a distance  increment,  the  electric  field  is 

j [ 

entirely  in  the  y direction  (defined  as  normal  to  the  waveguide  planes) , 

1 - 

at  the  output,  only  a Very  small  x and  z (the  axial  direction)  component 

? 4 

| ='  . 

has  been  generated.  In  our  experimental  system,  these  components  cannot 

55 
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be  measured.  However,  any  changes  in  the  y component  we  should  be  able  to 

1 I 

accurately  measure. 

The  TEM  waveguide  plates  are  sufficiently  separated  so  that 

avalanching  may  occur. 
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A Fourier  analysis  in  time  and  space  of  the  real  time  electric 

1 

field  would  reveal  that  this  electric  field  consists  of  short  and  long 
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wavelength  components.  The  long  wavelength  components  are  electromagnetic 

in  character,  and  will  propagate  in  the  axial  direction.  The  electric 

i 
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field  will  be  directed,  in  the  absence  of  a magnetic  field,  in  the  same 

IHj 

direction  as  the  electric  field  input  pulse,  namely,  normal  to  the  wave- 

n  • - 

guide  plates.  Thus,  the  plates  do  not  influence  the  propagation  of  these 

1111 

components . 

iSf- 

It  is  possible  that  there  will  be  short  wavelength  components. 

\mm 

not  predicted  by  the  AFWL  theory,  which  only  exist  in  the  nonlinear  regime. 

fillip 

These  are  electrostatic  in  character  (e.g.,  Bernstein  modes) ^ The 

14 


wavelength  of  these  components  is  typically  much  shorter  than  the  separation 
of  the  plates,  so  that  again,  the  plates  do  not  particularly  influence 
these  waves.  For  example,  for  Bernstein  modes,  collisional  damping  requires 
that  X,  the  wavelength  for  a Bernstein  mode, satisfy 


where  V is  defined  by  (12a)  and  f„  is  w / 2tt  and  the  cyclotron  frequency 
t c c ^ 

is  given  by  (12c).  For  our  system,  V is  approximately  5 x 10^  cm/sec  and 
for  the  scaling  that  we  have  anticipated  (S  = 310) , f is  300  MHz  if  we 
take  into  account  the  earth’s  magnetic  field.  This  leads  to 


X < .333  cm 


It  is  these  short  wavelength  waves  which  are  the  result  of  plasma  turbulence. 
Thus,  the  plates  do  not  change  the  quantitative  correspondence  between  the 


laboratory  model  and  the  ionosphere. 


IV.  General  Description  of  System  Components 

The  details  of  the  system  we  have  built  will  be  described  in  the 
next  section.  Here  we  give  a general  overview  of  this  system. 

The  system  should  contain  two  types  of  plasmas . The  one  which 
has  been  constructed  at  the  present  time  we  shall  refer  to  as  an  electron 
beam-produced  plasma  system,  while  the  other  we  shall  call  a discharge 
plasma  system.  These  plasmas  will  be  discussed  below,  as  well  as  the  TEM 
waveguide,  the  pulse  generator,  and  the  diagnostic  systems. 

A.  Electron  Beam-Produced  Plasma 

This  plasma  has  been  described  in  the  literature^’ , and  has 
been  recently  constructed  at  two  major  laboratories^^"’  In  Fig.  7 
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there  is  shown  a grid-cathode  arrangement  followed  by  a long  drift  space  at 
the  end  of  which  is  a collector.  The  grid,  which  is  close  to  the  cathode, 
accelerates  primary  electrons  emitted  by  the  cathode.  These  pass  through 
the  grid  with  energies  of  typically  50  ev. 

Now  the  entire  volume  of  the  system  contains  a gas  of  nitrogen  at 
a low  pressure.  Further,  an  external  solenoid  produces  an  axial  magnetic 
field.  The  primary  electrons  passing  through  the  grid-cathode  space  are  con- 
strained to  drift  along  a magnetic  field  line.  In  the  drift  space,  they 
ionize  gas  atoms,  producing  secondary  electrons  with  relatively  low  energies, 
typically  1 ev.  A secondary  electron  and  the  positive  ion  produced  when  the 
secondary  electron  is  removed  from  a neutral  atom,  also  move  along  magnetic 
field  lines.  Typically,  a primary  50  ev  electron  will  produce  several 
secondaries.  The  secondaries  have  approximately  a Maxwellian  velocity  dis- 
tribution. We  have  observed  that  if  the  collector  is  at  grid  potential,  the 
primary  electrons  make  one  pass  through  the  system  and  then  are  collected  by 
the  collector;  if  .he  collector  is  negative  with  respect  to  the  grid,  the 
primary  electrons  will  make  several  passes . 

Because  the  electrons  are  constrained  to  follow  magnetic  field 
lines,  the  plasma  cross-section  will  be  an  image  of  the  emitting  surface  of 
the  cathode.  Our  system  contains  a 3-inch  diameter  dispenser-type  cathode 
which  gives  a plasma  cross-section  whose  electron  density  has  cylindrical 
symmetry.  The  plasma  has  been  observed  to  nicely  fill  the  space  between 
the  TEM  waveguide  plates. 

We  shall  now  discuss  the  parameters  attainable  in  this  sytem.  A 
series  of  external  solenoids  produce  a magnetic  field  coaxial  with  the  axis 
of  the  plasma.  The  magnetic  field  has  a maximum  steady-state  value  of  300 
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gauss-  Within  the  plasma,  the  magnetic  field  is  constant  in  time  to  + 0.1%. 
The  minimum  magnetic  field  that  can  be  attained  while  still  confining  the 
plasma  is  about  15  gauss.  Thus,  for  a field  variation  of  15  to  300  gauss. 


the  variation  of  the  electron  cyclotron  frequency  uc/2tt  is  40  MHz  to  2.8C0 


GHz . I 

: 

The  electron  and  ion  number  density  depends  on  the  gas  pressure  ■ 

11  3 

and  on  the  grid-cathode  current . The  maximum  number  density  is  10  /cm  and  ? 

the  minimum  is  determined  by  diffusion  properties  of  the  electrons,  and  is 


6 3 

estimated  to  be  about  10  /cm  . This  gives  a range  for  the  electron  plasma 
frequency  u / 2tr  of  8.98  GHz  down  to  about  9 MHz. 

The  electron  momentum  transfer  collision  frequency  v depends  on  the 
kind  of  gas,  electron  temperature,  and  gas  pressure.  Nitrogen  can  be  used 
without  damaging  the  cathode.  We  can  generally  operate  at  pressures  in  the 
range  0.1  micron  to  20  microns.  For  He,  this  corresponds  to  a collision 
frequency  of  v of  10  MHz  to  100  KHz,  and  undoubtedly  higher  collision  fre- 
quencies than  10  MHz  would  be  possible.  For  it  is  roughly  the  same 
providing  that  the  electron  temperature  TeQ  is  not  a critical  value  of  2.2 
ev.  Then  the  collision  frequency  would  be  about  a factor  of  five  higher 


than  those  stated  above.  The  electron  temperature  T will  be  about  1 ev. 


and  is  not  under  our  control. 


B.  Discharge  Plasma 

This  plasma  is  produced  in  a large  pyrex  tube,  which  is  rectangular 
in  cross  section  where  the  TEM  waveguide  is  located  and  which  tapers  to  a 
round  cross  section  beyond  this  region. 

The  tube  is  open  at  both  ends.  A cathode  (1  inch  in  diameter. 


which  we  have  on  hand)  protrudes  into  the  round  end.  At  the  rectangular  end 


is  a collector.  By  grounding  the  collector  and  making  the  cathode  at  a 
negative  potential,  it  is  well  known  that  an  arc  discharge  can  be  established 
which  is  contained  by  the  glass  tube.  We  have  already  mocked  up  this  system 
using  a round  cross  section  tube  with  nitrogen  gas  and  found  that  an  arc 
discharge  can  easily  be  obtained. 

The  TEM  waveguide  plates  set  on  the  flat  surfaces  of  the  rectangu- 
lar portions  of  the  pyrex  tube.  They  are  external  to  the  discharge.  The 
details  of  construction  may  be  seen  in  Section  V. 

The  parameters  obtainable  in  this  type  of  plasma  correspond  to 
higher  pressures,  of  the  order  of  10  microns  or  higher  and  number  densities 
in  the  range  10^  to  10"^.  The  electron  temperature  is  about  the  same  as  in 
the  electron  beam-produced  plasma,  i.e.,  about  1 ev. 

C.  Pulse  Generator 

The  pulse  generator  consists  of  a mercury  reed  switch  in  series 
with  a DC  power  supply  of  voltage  at  the  desired  amplitude  of  the  EMP 
and  a coax  line  connected  to  the  TEM  waveguide  input.  It  is  opened  by  a 
magnetic  field  coil.  The  details  are  given  in  Section  V.  We  have  built  a 
device  which  has  a rise  time  of  100  picoseconds  and  a fall  time  which  we  can 
vary  by  the  size  of  a metal  slug  in  the  system.  Currently,  the  fall  time  is 
about  3 nanoseconds,  but  this  can  be  altered.  The  largest  amplitude  we  have 
introduced  into  the  plasma  loaded  waveguide  to  date  is  2000  volts.  We  have 
the  capability  of  increasing  this  to  nearly  10,000  volts.  We  do  not  know 
yet  if  arcing  will  occur  at  these  high  voltages. 

D.  TEM  Waveguide 

This  waveguide  consists  of  two  parallel  plates . Each  is  7 centi- 


meters wide,  and  they  are  separated  by  3 centimeters  wide  such  that  the 
characteristic  impedance  in  free  space  is  109  ohms.  In  the  electron 


beam-produced  plasma,  the  plates  are  1/4  inch  mesh  stainless  steel  screen 
while  in  the  discharge  plasma  they  are  strips  of  0.050  inch  thick  stainless 
sheet  which  will  be  screwed  to  the  top  and  bottom  of  the  pyrex  glass  dis- 


charge tube. 

For  either  type  of  plasma,  conical  tapers  provide  a broad-band 
transition  to  a twin-lead  transmission  line  which  passes  through  the  vacuum 
wall.  The  characteristic  impedance  of  109  ohms  is  maintained  through  the 
vacuum  wall  at  both  the  input  and  output. 

Outside  the  vacuum,  broad-band  tapers  provide  a transition  from 
the  109  ohm  twin-lead  line  to  a 50  ohm  coaxial  cable.  This  cable  attaches 
to  the  pulse  generator  at  the  input  and,  through  a suitable  termination,  to 
the  sampling  scope  at  the  output. 

E.  Electrostatic  Analvzer  Probe 


In  the  presence  of  an  IMP,  the  electrons  will  be  heated  such  that 

their  average  energy  far  exceeds  1 ev.  We  are  expecting  peaks  of  average 

energy  exceeding  100  ev.  Furthermora,  this  tvpe  of  heating  may  last  only  a 

short  time;  after  the  electrons  are  heated  by  the  EMP,  and  the  pulse  has 

propagated  out  of  the  waveguide,  they  will  lose  their  energy  through 

collisions  with  neutrals  or  will  diffuse  out  of  the  system.  The  plasma  will 

then  return  to  equilibrium  T . The  electrons  will  thus  be  at  a high  energy 

for  a short  time,  probably  for  considerably  less  than  1 microsecond. 

In  order  to  measure  the  energy  of  these  hot  electrons,  we  have 

built  an  electrostatic  analyzing  probe.  This  probe  has  been  developed  by 
[13] 

Mix,  et  al.  , who  used  this  device  for  measuring  the  distribution  ot 
electron  energies  in  an  electron  beam.  The  application  of  their  device  to 
measuring  electron  energies  in  a plasma  at  energies  up  to  300  ev  has  been 
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developed  by  M.  IWKolaD,  et  al,  In  particular,  a student  of  M, 

Porkolab  has  adapted  this  device  to  a machine  very  similar  to  our  machine 
at  The  University  of  Arizona. 

This  device  not  only  determines  weak  electron  energy,  but  will 
also  permit  a determination  of  the  energy  distribution  function.  Measure- 
ment of  the  energy  distribution  function  is  performed  on  a time  scale  of 
about  5 nanoseconds.  From  this,  number  density  and  drift  velocity  should 
1 e determinable. 

F . Electromagnetic  Cavity  and  Langmuir  Probes 

A standard  electromagnetic  cavity  is  being  used  to  measure  ambient 
[14] 

number  density  . To  measure  ambient  electron  temperature  Teo,  we  are 
using  two  Langmuir  probes,  one  which  can  be  moved  radially  into  the  plasma, 
and  the  other  which  can  be  moved  longitudinally  along  the  TEM  waveguide. 

The  Langmuir  probe  is  pulsed.  A voltage-current  plot  on  an  oscilloscope  is 
obtained  by  applying  a sawtooth  voltage  to  the  probe  where  the  voltage  is 
swept  from  zero  to  maximum  in  3 ys.  As  is  well  known,  from  the  slope  of 
this  voltage-current  plot,  TeQ  may  be  determined. 

V . Description  of  the  Equipment  That  Has  Been  Built 
Under  Air  Force  Contract 

A.  General  Layout  of  System 

A schematic  drawing  of  the  electron-beam  produced  plasma  system  is 
shown  in  Fig.  8a.  The  horizontal  portion  of  the  vacuum  system  is  about  7 
feet  long.  A photograph  of  the  completed  machine  is  given  in  Fig.  8b. 

B.  Components 

1.  Vacuum  System  and  External  Vacuum  Envelope 

The  vacuum  envelope  is  of  304  stainless  steel,  about  8 inches  0D 
and  approximately  6 feet  long.  There  is  a 45°  side  port  through  which 
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emerges  the  output  of  the  TEM  waveguide.  The  vacuum  system  is  a conventional 
diffusion  pump  backed  by  a fore  pump.  Pressures  at  the  remote  end  of  the 
chamber  of  order  of  10  ^ torr  have  been  achieved. 

2.  Pressure  Measuring  System 

We  use  a thermocouple  gauge  to  measure  air  pressure  down  to  1 

..7 

micron,  and  ion  gauges  for  air  pressure  down  to  10  torr.  These  gauges  are 
used  in  evacuating  the  system  initially.  Nitrogen  is  admitted  through  a 
needle  valve,  and  its  pressure  in  the  system  is  monitored  by  Shultz-Phelps 

gauges  which  can  measure  nitrogen  pressure  in  the  range  10  ^ to  1 torr. 

♦ 

3.  Cathode,  Grid,  and  Collector 

The  cathode  is  a dispenser-type  cathode  approximately  3 inches 
in  diameter.  A grid  of  close  spaced  tungsten  wires  is  placed  about  3/^  inch 
in  front  of  the  cathode.  Emission  of  electrons  from  the  cathode  surface  is 
controlled  by  the  voltage  difference  between  the  grid  and  cathode.  The 
resulting  beam  of  electrons  interacting  with  the  nitrogen  gas  produces  the 
plasma.  The  ambient  number  density  of  the  plasma  is  controlled  by  the 
electron  beam  current  which  in  turn  is  controlled  by  the  grid-cathode  volt- 
age. We  have  found  that  the  most  stable  operation  can  be  obtained  with  the 
grid  somewhat  above  ground.  The  waveguide  and  input  conical  transition  serve 
as  the  collector.  These  elements  are  grounded. 

4.  Waveguide  Assembly 

The  TEM  waveguide  parallel  plates  are  made  from  stainless  steel 
1/4  mesh  screen  welded  onto  a supporting  track.  Each  plate  is  7 centimeters 
wide,  and  they  are  separated  by  3 centimeters.  The  waveguide  is  about  1 
meter  long  with  the  axis  at  the  output  at  45°  from  the  axis  at  the  input  as 
shown  in  Fig.  8.  We  have  found  that  the  45°  bend  introduces  negligible 
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reflections.  At  each  end  of  the  waveguide  is  a conical  transition  . These 
transition's  are  identical.  Each  is  8 inches  long.  A transition  transforms 
geometrically  from  the  TEH  parallel  plate  waveguide  to  a twin-lead  trans- 
mission line.  The  twin  lead-transmission  line  is  enclosed  in  a circular 
cylindrical  stainless  tube  and  embedded  in  Torr-seal^1^  (a  vacuum  sealing 
dielectric  with  a relative  permittivity  of  4.0).  The  twin-lead  transmission 
thus  passes  through  the  vacuum  wall.  The  characteristic  impedance  of  all 
portions  of  the  transmission  line  is  109  ohms^^.  At  both  the  input  and 
output % identical  tapers  transform  the  109  ohm  twin-lead  to  50  ohm  coax. 

A cross-section  of  a taper  is  shown  in  Fig.  9. 

5 . Pulse  Generator 

To  generate  our  fast  rise-time,  high  voltage  pu]ses,  we  are  using 
a modified  charged-line  pulse  generator,  the  construction  of  which  is  shown 


in  Fig.  10.  Our  pulse  can  be  characterized  by  three  parameters:  (1)  rise- 

time, (2)  fall-time,  and  (3)  peak  value.  The  rise-time  we  would  like  to 
minimize.  The  fall-time  we  would  like  to  be  able  to  vary.  The  rise-time 


would  be  minimized  when  we  have  minimized  all  series  inductances.  In  em- 
bedding the  switch  inside  the  inner  conductor,  we  have  minimized  the 
inductance. 

Maximizing  the  peak  value  requires  that  the  outer  radius  of  the 
inner  conductor  approach  the  inner  radius  of  the  outer  conductor  since 


V , , . = V.  Z./Z.  + Z 

o(peak)  mil  o 


where 


Z^  = 50  ohms 

Z = 42.5  £n  b/a 
o 


, *--•  - n t 


b = inner  radius  of  outer  conductor 


a = outer  radius  of  inner  conductor 


If  we  were  to  do  this,  the  capacitance  per  centimeter  would 


approach  infinity. 


C = 1.11  to  b/a  / picofarad/c 


The  fall-time  is  directly  related  to  the  length  and  diameter  of 


the  inner  conductor: 


Fall-time  = 2.2  • R » C 


R = Z_  + Z 
1 o 


C = 1.11  in  b/a  f picofarad/* 


The  minimum  length  of  the  inner  conductor  that  we  can  use  is 
approximately  1 inch,  the  length  of  the  reed  switch  plus  lead  length. 

With  all  the  above  in  mind,  a 5 nanosecond  fall-time  pulse  genera- 
tor was  built.  This  gives  a maximum  output  voltage  of  87%  of  the  input 
voltage. 


'I 


The  measured  results  were  as  follows: 


Rise-time  160  picoseconds 


Fall-time  = 5 nanoseconds 


V „ = 80%  V. 

out  in 


The  magnetic  field  used  to  operate  the  reed  switch  was  obtained 
by  winding  approximately  150  turns  of  wire  around  the  pulse  generator  and 


then  applying  a voltage  to  the  coil.  The  maximum  operating  frequency  of 
the  reed  switch  is  about  200  Hz. 
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6 . Magnets 


Seven  magnets  surround  the  system  as  may  be  seen  in  Fig.  8.  In 
addition,  a mirror  magnet  surrounds  the  input  taper.  The  function  of  this 


latter  magnet  is  to  reflect  electrons  which  enter  the  V-shaped  region  of  the 
conical  transition.  Since  an  electron  travelling  toward  the  cathode  will  be 
reflected  in  the  region  between  the  grid  and  cathode,  electrons  will  tend  to 
be  trapped  longitudinally,  i.e.,  they  will  be  reflected  at  the  two  ends. 

The  current  in  each  of  the  eight  magnets  is  adjustable  by  means  of 
potentiometers.  We  have  a computer  code  which  tells  us  how  to  set  the  current 
in  each  magnet  for  any  desired  magnetic  field  variation  on  the  axis  of  the 


M 


vacuum  chamber. 


C.  Discharge  Plasma  System 

(Not  constructed  as  of  Jan.  15,  1975) 

The  discharge  plasma  system  is  shown  schematically  in  Fig.  11.  A 
glass  tube,  shown  in  Fig.  12,  rectangular  in  cross  section  at  one  end  and 
circular  at  the  other  is  mounted  in  the  vacuum  chamber.  In  place  of  the 
screen  TEM  waveguide,  two  plates,  made  from  0.050  inch  thick  stainless  sheet, 
with  identical  dimensions  to  the  screen  TEM  waveguide  (7  cm.  wide,  3 cm.  apart, 
and  the  same  longitudinal  dimensions)  are  screwed  to  the  top  and  bottom  ex- 
ternal surfaces  of  the  gJass  pyrex  tube.  All  other  parts  of  the  transmission 
line  system  are  interchangeable  with  the  screen  TEM  waveguide,  e.g.,  conical 
transitions,  twin-lead  line,  and  external  tapers  to  coax. 
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The  3 inch  cathode  is  replaced  by  a smaller  1 inch  cathode  which 
protrudes  into  the  round  portion  of  the  glass  tube.  This  may  be  seen  in 


Fig.  11. 


An  arc  discharge  can  be  obtained  between  the  cathode  and  the 


input  conical  transistion  which  now  acts  as  an  anode  for  the  discharge.  We 
have  mocked  up  this  system,  and  have  verified  that  it  will  work. 


D. 


Diagnostics 

1.  EMP  Propagation  Analyzer 
% 

To  determine  what  effect  our  plasma  will  have  on  the  EMP,  we  will 
use  the  system  shown  in  Fig.  13.  The  switching  circuit  produces  a high 
voltage  fast  rise- time  pulse  with  an  exponential  decay  as  discussed  above. 

The  pulse  is  then  split  by  a power  divider  enabling  us  to  view  the  signal 
before  it  has  entered  the  plasma  system.  Since  sampling  scopes  typically 
take  120  nanoseconds  to  respond  to  a trigger  signal,  we  required  a 120 
nanosecond  delay  line  between  the  pulse  and  the  vertical  input  of  the 
oscilloscope.  The  pulse  propagating  through  the  plasma  will  also  require 
a delay  line  of  120  nanoseconds  before  going  into  the  vertical  input  of  the 
oscilloscope.  We  now  have  a trace  of  the  input  to  the  plasma  and  a trace 
of  the  pulse  after  it  has  been  propagated  by  the  plasma. 

2.  Electrostatic  Analyzer  Probe 

The  Electrostatic  Analyzer  Probe  has  as  its  function  to  measure 
the  electron  distribution  function.  It  will  certainly  allow  detection  of 
100  ev  electrons.  Under  our  present  design,  we  can  measure  the  distribu- 
tion function  either  at  the  time  the  pulse  passes  through  the  waveguide  or 
100  or  200  picoseconds  later.  The  probe  is  longitudinally  movable,  so  that 
the  distribution  function  as  a function  of  longitudinal  distance  can  be 
obtained,  if  desired.  It  can  be  used  with  either  plasma  system. 

Figure  14  is  a schematic  of  how  we  obtain  the  distribution  function 
(velocity  or  energy)  of  the  plasma  after  it  has  been  excited  by  the  pulse. 
Power  supply  A provides  an  electric  field  such  that  ions  are  excluded  from 
the  electrostatic  probe.  By  varying  power  supply  B,  we  obtain  the  integral 
of  the  electron  distribution  function.  This  is  because  supply  B sets  up  an 
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electric  field  in  the  analyzer  that  repels  electrons  with  energies  lower 


than  the  supply  voltage-  Therefore,  the  number  of  electrons  reaching  the 


collector  is  directly  proportional  to  the  supply  voltage.  The  audio 


oscillator  and  lock-in  amplifier  serve  as  a means  of  differentiating  the 


collector  current,  thereby  giving  us  the  distribution  function. 


3.  Electromagnetic  Cavity 


The  Electromagnetic  Cavity  allows  the  measurement  of  ambient 

i — it...  i _ .-i i [l^J  -t  - i„  j 


average  electron  number  density  in  the  plasma  column 


It  is  designed 


to  be  compatible  with  either  plasma  system.  In  the  magnetically  confined 


plasma,  it  is  longitudinally  movable  when  there  is  no  waveguide  in  the 


vacuum  system,  and  it  will  be  used  in  this  mode  to  initially  determine 


longitudinal  number  density  variations.  With  the  waveguide  in  place,  it 


remains  fixed  near  the  cathode  and  so  can  measure  electron  density  at  only 


that  longitudinal  position.  It  can  measure  ambient  number  densities  as  low 


7 3 

as  10  /cm  . In  order  to  determine  that  we  have  number  densities  in  the 
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range  of  10  /cm  , we  will  have  to  extrapolate  our  cathode  current  and  voltage 


7 3 

data  from  the  region  of  10  /cm  and  above. 


We  use  the  microwave  cavity  to  obtain  an  average  number  density 


of  the  plasma  column.  Basically,  we  measure  the  resonant  frequency  shift 


of  the  TMq-^q  mode  (lowest  cutoff  frequency)  when  a plasma  is  present.  This 


frequency  shift,  Af,  is  related  to  the  plasma  frequency  f by: 


f2  = M Aff  (1  + v2/f2) 
P 


where 


M is  a dimensionless  constant  which  is  a function  of  the 


ratio  b/a  (rig.  15) 


v is  the  collision  frequency 

f is  the  measured  resonant  frequency  of  the  cavity 


with  the  plasma  present. 


2 2 

For  our  experiment  v /f  « 1. 
Fig.  16. 


The  experimental  apparatus  is  shown  in 


4.  Langmuir  Probe  Diagnostics 

The  Langmuir  Probe  allows  absolute  measurement  of  ambient  electron 
temperature.  It  also  allows  relative  measurement  of  ambient  number  density. 
The  probe  is  pulsed  over  a 3 microsecond  interval,  so  "ambient"  refers  to 
that  interval.  Thus,  if  there  is  a slow  variation  of  either  number  density 
or  temperature  compared  to  3 microseconds,  the  Langmuir  probe  will  be  capable 
of  measuring  these  quantitites  as  a function  of  time.  The  distribution 
function  must  approximate  Maxwellian  (probably  unlikely  except  in  the  truly 
ambient  case)  for  the  measurements  to  have  meaning  that  we  can  interpret. 

Since  the  cavity  allows  absolute  measurement  of  number  density  at 
a fixed  point,  this  piece  of  information  can  be  used  to  convert  relative 
number  density  measurement  to  absolute  measurements.  The  error  in  doing  this 
is  unknown  at  the  present  time. 

We  anticipate  using  Langmuir  probes  to  obtain  electron  temperature 
and  relative  number  density.  There  will  be  two  probes,  as  shown  in  Fig.  17. 
One  probe  will  be  used  to  obtain  radial  measurement  [n(r),  T(r)],  and  can 
be  inserted  at  two  locations  into  the  plasma.  The  other  probe  will  be  used 
for  axial  measurements  over  the  entire  length 


[n(z),  T(z) ] 


of  the  structure.  These  probes  can  be  used  with  either  the  electron  beam- 
produced  or  discharge  plasma  systems.  Due  to  the  presence  of  the  magnetic 
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field,  it  will  be  difficult  tc  obtain  accurate  number  densities  from  these 
probes.  However,  relative  density  profiles  can  be  obtained,  and  compared 
with  cavity  measurements  to  obtain  absolute  number  densities.  Further,  the 
magnetic  field  anticipated  should  not  appreciably  affect  temperature 
measurements  using  the  probe. 

Figure  18  is  a schematic  of  our  Langmuir  probe  instrumentation. 

This  equipment  has  been  built  and  is  presently  operating.  We  have  measured 

9 11  3 

number  densities  of  10  to  10  elec/cm  . A larger  probe  will  be  used  for 

7 8 

diagnostics  at  lower  number  densities  (10  to  10  ) . 


VI.  Preliminary  Results 

In  this  section,  we  present  results  obtained  from  the  plasma 

machine,  as  of  January  15,  1975. 

A.  Cold  Testing  the  Waveguide  with  Time 
Domain  Reflectometer  (TDR) 

To  determine  qualitatively  whether  -he  transitions  will  be  of 
beneficial  help  to  pulse  propagation,  a TDR  analysis  was  run  on  the  system. 
The  setup  used  is  shown  in  Fig.  19a. 

If,  instead  of  the  transitions  we  had  a direct  mismatch  at  50  to 
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109  ohms,  we  would  expect  a TDR  plot  such  as  shown  in  Fig.  19b.  If  the 
tapers  provided  a perfect  match,  we  would  expect  a horizontal  line  TDR  plot. 
This  would  correspond  to  a reflection  coefficient  of  0.  Since  the  tapers 
have  a reflection  coefficient  that  varies  with  frequency,  being  higher  at 
low  frequencies  than  high  frequencies,  we  should  be  able  to  qualitatively 
predict  our  TDR  pattern. 

Since  more  low  frequencies  will  be  reflected  than  high  frequencies, 
we  would  expect  the  steps  on,  the  TDR  plot  to  be  gradual  rises  rather  than 


sharp  rises.  If  this  is  the  case,  then  our  tapers  will  have  served  their 
purpose.  This  statement  will  be  explained  later. 

The  TDR  plot  for  our  system  is  shown  in  Fig.  19c.  As  can  be  seen, 
there  are  small  sharp  reflections  along  with  large  gradual  reflections.  The 
small  reflections  have  reflection  coefficients  of  about  0.08  and  are  too 
small  to  worry  about.  The  large  gradual  reflections  are  as  predicted. 

The  purpose  behind  using  the  tapers  was  to  obtain  a higher  voltage 
through  the  waveguide  with  as  little  effect  on  the  rise-time  of  the  pulse 
as  possible.  If  the  tapers  hadn't  been  used,  we  could  only  have  expected 
approximately  60%  of  the  input  pulse  voltage  to  propagate  down  the  wave- 
guide, and  only  60%  of  that,  or  35%  of  the  input  pulse,  propagating  out  of 
the  system.  This  would  mean  that  we  would  have  approximately  24%  of  the 
input  pulse  being  reflected  back  down  the  waveguide,  which  is  not  a very 
desirable  situation.  Since  the  reflection  coefficient  of  the  tapers  is 
lower  at  high  frequencies,  and  since  we  are  using  high  frequency  pulses,  our 
pulse  peak  amplitude  will  be  greater  than  if  we  had  had  a straight  mismatch. 
The  "average"  reflection  coefficient  of  each  taper  can  be  calculated.  We 
have  found  that  the  taper  allows  between  70  and  80%  of  the  pulse  voltage  to 
propagate  in  the  waveguide. 

B.  Pulse  Characteristics  of  the  System  without  Plasma 

The  system  configuration  is  shown  in  Fig.  20.  The  pulse  generator 
output  is  divided,  half  traveling  through  the  system  while  the  other  half  is 
used  to  trigger  the  scope  and  observe  the  input  pulse. 

The  rise-times  of  the  different  components  are: 
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Rise-time 


Equipment 
Oscilloscope 
Delay  Lines 
Power  Dividers 
Attenuators 


100  picoseconds 
120  picoseconds 
25  picoseconds 
25  picoseconds 


Since  we  are  working  with  pulses  with  rise-times  close  to  that  of 
the  measuring  equipment,  we  must  correct  for  these  rise-times.  This  is  done 
by  the  following: 

2 — 2 i 

(measured  rise-time)  - £ (rise-time  of  equipment) 

The  measured  rise-time  of  the  input  and  output  pulses  was 
respectively  200  ps  and  450  ps.  Therefore,  the  actual  worst  case  rise-times 
of  the  input  and  output  pulses  are  respectively  160  ps  and  425  ps.  Assuming 
that  the  rise-time  of  the  waveguide  is  very  small  and  the  rise-times  of  the 
transitions  are  equal,  both  very  good  approximations,  we  calculated  the  rise- 
times  of  the  tapers  to  be  300  ps.  Knowing  this,  the  rise-time  of  the  pulse 
in  the  waveguide  is  calculated  to  be  about  340  ps. 

A plot  of  the  input  and  output  peak  voltages  as  a function  of  the 
voltage  across  the  switch  is  shown  in  Fig.  21.  As  can  be  seen,  the  output 
voltage  is  down  about  30%.  This  is  due  to  the  rise-time  of  the  tapers  being 
slower  than  the  rise-time  of  the  pulse.  The  maximum  output  voltage  of  the 
system  shown  in  Fig.  21  is  about  700  volts.  This  value  could  be  doubled  by 
removing  the  power  divider  at  the  input  and  connecting  the  pulser  directly 
to  the  taper.  The  trigger  signal  could  then  be  obtained  from  the  output 
pulse.  The  pulse  shape  was  found  to  be  independent  of  the  voltage  as 
would  be  expected. 
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The  following  is  a summary  of  the  range  of  values  we  are  presently 
a'oie  to  obtain: 


Input  rise-times 
Output  rise-time 
Input  fall-time 
Output  fall-time 
Input  voltage 
Voltage  in  waveguide 
Output  voltage 


160  picoseconds 
450  picoseconds 
variable 

same  as  input  fall-time 
0 - 2000  volts 

(0  -1600)  volts  (approximate) 
0 - 1400  volts 


C.  Preliminary  Study  of  Pulse  Propagation 
Through  Plasma  Filled  Waveguide 

We  have  some  very  preliminary  results  from  propagating  an  EMP  type 

pulse  along  the  TEM  waveguide  filled  with  plasma.  The  only  diagnostic  which 

was  installed  was  the  radial  Langmuir  probe.  This  allowed  only  a crude 

measurement  of  ambient  number  density  n , which  we  estimated  to  be  about 

eo 

9 3 18  3 

10  /cm  . The  ambient  neutral  density  was  5.16  x 10  /cm  corresponding  to  a 

-3 

pressure  of  1.6  x 10  torr.  In  the  absence  of  plasma,  the  output  pulse 

appeared  as  shown  in  Fig.  22  a,  c,  and  e.  The  peak  amplitude  of  the  output 

pulse  is  700,  540,  and  260  volts,  respectively.  These  pulses  should  be  very 

nearly  the  pulses  that  appear  at  the  input  to  the  waveguide. 

Regarding  Fig.  22a  in  an  input  pulse  to  the  waveguide,  when  a 

plasma  is  present,  the  corresponding  output  appears  as  in  b.  Similarly,  for 

input  pulses  c and  e,  the  output  pulses  are  d and  f,  respectively. 

In  the  output,  initially,  a ringing  effect  occurs  with  a period 

of  1.4  nanoseconds  (ns)  or  a frequency  of  710  MHz.  If  one  assumes  that  this 

frequency  is  the  plasma  frequency,  then  ngo  that  is  predicted  by  this  fre- 
9 3 

quency  is  6.3  x 10  /cm  , which  is  in  rough  agreement  witn  the  Langmuir  probe 
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estimate  of  n . 
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During  these  measurements,  the  confining  magnetic  field  was  low 
enough  that  a further  reduction  of  the  magnetic  field  did  not  alter  the 
results. 

After  about  1-1/2  periods  of  ringing,  the  output  pulse  amplitude 
increased  for  about  5 ns  and  thereafter  went  slowly  to  zero  after  about  an 
additional  10  ns.  The  reason  for  this  behavior  is  unknown. 

VII.  Recommendations  for  Future  Work 

By  the  end  of  the  technical  phase  of  this  work  (January  15,  1975), 
the  machine  itself  had  been  completed.  The  diagnostic  subsystems  had  also 
been  completed,  with  the  exception  of  the  electrostatic  analyzer  probe,  but 
none,  except  the  radial  Langmuir  probe,  had  been  introduced  into  the  machine. 
Thus,  the  problems  of  verifying  that  these  subsystems  will  do  their  job  and 
work  harmoniously  with  each  other  still  remain  to  be  worked  out.  In  addition, 
some  preliminary  data  have  been  taken,  although  certainly  not  enough  to  shed 
much  light  on  the  validity  of  the  AFWL  theory. 

We  present  here  our  recommendations  for  the  directions  that  this 
project  should  take,  should  continued  funding  become  available.  Specifically, 
we  would  recommend  the  following: 

1.  We  should  introduce  the  various  diagnostic  subsystems  (electro- 
static analyzer  probe,  cavity,  Langmuir  probes)  into  the  plasma  machine,  and 
test  each  subsystem.  Each  must  work  in  concert  with  the  rest.  There  will 
be  some  time  required  to  optimize  the  entire  system.  For  example,  at  the 
present  time,  we  have  poor  control  over  nitrogen  gas  pressure  because  of  the 
high  pumping  speed  of  the  diffusion  pump.  We  anticipate  that  this  can  be 
overcome  by  a baffle  that  can  be  placed  near  the  cathode. 
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2.  Construct  the  discharge  plasna.  The  only  component  nissing  at 
the  present  time  is  the  pyrex  tube  shown  in  Fig-  12-  We  have  a bid  on  this 
item  indicating  that  its  fabrication  will  cost  approximately  $500-  All 
other  parts  are  on  hand.  This  system  has  been  designed  to  be  compatible 
with  all  diagnostic  subsystems-  Some  optimization  will  be  required - 

3-  The  complete  plasma  machine  should  then  be  used  to  obtain  extensive 
experimental  aata  or  the  effect  of  plasma  on  the  propagation  of  an  EMP 
through  a plasma,  as  well  as  the  effect  of  the  pulse  on  the  plasna. 

Specifically,  the  following  quantities  should  be  measured  as  a 

function  of  the  independent  variables  listed  below.  These  quantities  are: 

maximum  energy  of  electrons.  Te  which  should  be  measured  as  a function 

max  j 

of  longitudinal  position;  distribution  function,  f(T  ),  to  be  measured  at 

e 

discrete  longitudinal  points;  output  electric  field  from  the  TEH 
transmission  line,  E(t),  to  be  measured  in  real  time  using  a sampling 
scope-  The  independent  variables  are:  a)  input  EMP  pulse  maximum  ampli- 

tude, Eq,  b)  EMP  input  pulse  rise-time,  t^.:  c)  the  input  pulse  fall-time, 

t,;  d)  ambient  number  densitv,  . ; e)  electron  neutral  collision  fre- 

i eo 

quency  v.  The  independent  parameters  should  be  chosen  such  that  an 

incremental  distance  in  the  D-region  can  be  simulated. 

4.  We  should  make  a careful  measurement  of  T^o  with  the  Langmuir 

probe.  It  is  expected  that  this  will  be  close  to  1 ev.  Using  this  value 

of  T , we  need  to  have  a comouter  run  made  with  T set  equal  to  this  value, 
eo  • eo 

The  run  should  be  made  for  one  distance  increment,  and  u,E,t,tf,  and  v 

p oT  r f 

chosen  such  that  the  scaled  values  are  attainable  in  our  laboratory  system 
and  the  theoretically  predicted  results  show  electron  temperatures  approach- 
ing 100  ev.  The  results  from  this  run  can  easily  ba  scaled  to  predict  the 


experimental  results  for  a range  of  in  our  lal  oratory  plasma,  where  the 
scaling  is  to  be  done  according  to  (27). 

5.  In  order  to  understand  the  discrepancies  between  theoretical  and 
experimental  E(t),  a computer  code  should  be  obtained  which  will  Fourier 
transform  E(t)  into  the  frequency  domain  E(io).  Thus,  for  example,  if 
discrete  frequency  components  are  dissimilar,  this  dissimilarity  will  be 
quickly  and  quantitatively  perceived. 

6.  Eq  may  have  a significant  effect  on  f(w)  and  Temax.  This  should 
be  carefully  investigated. 

7.  The  experimental  system  should  be  modified  to  include  probes  which 

will  perform  the  following  functions:  a)  measure  maximum  number  density  due 

to  avalanching;  b)  measure  average  electric  field  as  a function  of  longitudi- 
nal distance  (the  average  electric  field,  according  to  the  theory,  should 
change  nearly  linearly  with  distance);  c)  measure  the  presence  of  electro- 
static waves  (these  are  not  predicted  by  the  theory  but  may  be  excited  by 


the  pulse) . 

8.  Summarize  the  correspondence  between  theory  and  experimental 
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Figure  1.  Geometry  of  one  dimensional  calculation  of  pulsed 
electromagnetic  fields  through  the  D-region. 


Figure  7.  Schematic  diagram  of  the  beam  plasma  machine  operation. 


Figure  8A.  A schematic  of  the  electron-beam  produced  plasma  systei 


Figure  9.  Impedence  tapes  which  mutches  109  ohm  twin-lead  to 
50  ohm  coax. 


Figure  10.  The  modified  charge-line  puJse  generator  used  to  generate 
the  fast  rise-time,  high  voltage  pulses. 


Schematic  of  the  EMP  propagation  analyzer. 


Figure  15.  The  electrostatic  cavity  used  for  the  determination  of 
electron  nimber  density. 


Figure  16.  Schematic  of  the  electrostatic  cavity  system  used  for 
the  determination  of  electron  number  density. 


(o)  IrfPUT  PULSE;  PULSE  PEAK  700  VOLTS.  VERTICAL  SCALE.  200  V/cm. 


(c> 


INPIF  PULSE;  PULSE  PEAR,  540  VOLTS. 


VERTICAL  SCALE,  200  V/ cm. 


(d)  OUTPUT  PULSE  CORRESPOND. -IS  TO  (c).  VERTICAL  SCALE,  100  V/cm. 


Figure  22. 


(continued) 
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ATTN:  L.  Loughlin,  Library 

Martin  Marietta  Aerospace 

ATTN:  End.  Lib.,  M.  C.  Griffith,  MP-30 

Martin  Marietta  Corp. 

ATTN:  D.  M.  Newell,  MS-S0454 

Maxwell  Lab. 

ATTN:  Dr.  V.  Fargo 

McDonnell  Douglas  Corp. 

ATTN:  C.  G.  Polak,  E-411-107-1-D2 

McDonnell  Douglas  Corp. 

ATTN:  R.  J.  DeBattista 
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McDonnell  Douglas  Corp. 

ATTN:  Library 

McDonnell  Douglas  Corp. 

ATTN:  T.  J.  Lundregan 

Mission  Research  Corp. 

ATTN:  D.  E.  Merewether 

Mission  Research  Corp. 

ATTN:  Dr.  C.  L.  Longmire 

Mitre  Corp. 

ATTN:  Library 

Mitre  Corp. 

ATTN:  L.  D.  Bickmoke 
Motorola,  Inc. 

ATTN:  TIC,  A.  J.  Kordalewski 

Northrop  Corp. 

ATTN:  J.  M.  Reynolds 

Northrop  R&T  Center 
ATTN:  Library 

Northrop  Corp. 

ATTN:  B.  T.  Ahlport 

Pac-Sierra  Research  Corp. 

ATTN:  Sec.  Ofc. 

Perkln-Elmer  Corp. 

ATTN:  T.  D.  Jones 

Philco-Ford 

WDLD 

ATTN:  Library 

Philco-Ford  A&C  Ops. 

ATTN:  Dr.  L.  H.  Linder 

Phys.  Internat. 

ATTN:  Dr.  F.  Ford 

Pulsar  Assoc. 

ATTN:  C.  Jones 

R &•  D Assoc. 

ATTN:  Dr.  W.  J.  Karzas 
ATTN:  R.  R.  Schaefer 

Rand  Corp. 

ATTN:  Dr.  C Crain 

Raytheon  Co. 

ATTN.  G.  Joshi,  Rad.  Sys.  Lab. 

Raytheon  Co. 

' ATTN:  J.  R.  Wcckback 
RCA  Corp. 

ATTN:  E.  Van  Keuren,  13-5-2 

RCA  Corp. 

ATTN:  E.  Daly 

RCA  Corp. 

ATTN:  Dr.  G.  Bruckor 
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Rockwell  Int.  Corp. 

ATTN:  TIC,  BA08 

Rockwell  Int.  Corp. 

ATTN:  J.  F.  Sexton,  CA31 

Sanders  Assoc. 

ATTN:  R.  G.  Despathy,  Sr.,  PE,  1-6270 

Science  Applications,  Inc. 

ATTN:  Dr.  F.  Tesche 

Science  Applications,  Inc. 

ATTN:  N.  R.  Byrn 

Science  Applications,  Inc. 

ATTN:  Dr.  R.  Beyster 
ATTN:  Dr.  G.  Stuart 

Science  Applications,  Inc. 

ATTN:  Dr.  B.  Goplen 

Simulation  Physics,  Inc. 

ATTN:  R.  G.  Little 
ATTN:  W.  A.  Sculler 

Simulation  Physics,  Inc. 

ATTN:  A.  J.  Armini 

The  Singer  Co . 

ATTN:  Dr.  A.  A.  Wittelcs,  M.S.  12-5820 

Sperry  Rand  Corp. 

UNIVAC  Div. 

ATTN:  J.  A.  Inda,  M.S.  5451 

Sperry  Rand  Corp. 

ATTN:  C.  Craig,  EV 

Stanford  Research  Institute 

ATTN:  Dr.  R.  A.  Armistead 
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Stanford  Research  Institute 
ATTN:  H.  Carey 

Sundstrand  Corp. 

ATTN:  C.  White,  Dept.  763SW 

Systems,  Science  &•  Software 
ATTN:  A.  F.  Klein 

Systron-Donner  Corp. 

ATTN:  H.  D.  Morris 

Texas  Instruments,  Inc. 

ATTN:  Dr.  D.  J.  Manus,  M.S.  72 

TRW  Systems  Group 

ATTN:  J.  E.  Dahnke 

TRW  Systems  Group 

ATTN:  W.  H.  Robinette 
ATTN:  TIC,  S-1930 

UAC 

Norden  Div. 

ATTN:  C.  Corda 

University  of  Denver 
Colorado  Seminary 

ATTN:  R.  W.  Buchanan 

University  of  Arizona 
Engineering  Research  Laboratories 
5 cy  ATTN:  Dr.  Robert  Carlile 

Varian  Assoc. 

ATTN:  Dr.  H.  R.  Jory,  A-109 

Weslinghouse  Electric  Corp. 

ATTN:  H.  P.  Kalapaca,  M.S.  3525 
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